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Abstract 



CN 

\q \ We present the first concrete evidence for the classical stability of 

■ vortons, circular cosmic string loops stabilized by the angular momen- 

tum of the charge and current trapped on the string. We begin by 
\ summarizing what is known about vorton solutions and, in particular, 

CO ■ their analytic stability with respect to a range of radial and nonradial 

perturbations. We then discuss numerical results of vorton simulations 
^ \ in a full 3D field theory, that is, Witten's original bosonic supercon- 

i ducting string model with a modified potential term. For specific pa- 

rameter values, these simulations demonstrate the long-term stability 
<-j ! of sufficiently large vorton solutions created with an elliptical initial 

ansatz. 

X ; PACS: 98.80.-k;74.60.Jg 

1 Introduction 



Topological defects are known to play an important role in many physical 
contexts and they may also impact cosmology, the relativistic setting for this 
present study (for a review, see ref. [I]). Amongst the possible cosmic defects, 
vortex-strings have a prominent position because they naturally evolve into 
a scale-invariant configuration, therefore avoiding analogues of the monopole 
problem of cosmological domination. These strings might be responsible 
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for a variety of astrophysical phenomena, such as cosmic rays, gravitational 
wave radiation, or gravitational lensing. The richness of their phenomenology 
comes in part from the possibility of additional internal structure, making 
them superconducting ref. [2]. 

These superconducting strings have been widely studied, and this article 
follows a companion paper [Hj, which investigated in detail the model studied 
here. Previously, we considered analytic criteria for current and charge sta- 
bility on these strings, confirming this with numerical simulations. This work 
allowed us to predict the equilibrium state of a cosmic string loop stabilized 
by the angular momentum of this charge and current, while also setting new 
limits on the available range of parameters. Let's mention that the potential 
existence of vorton solutions has been postulated recently in other important 
physical contexts such as QCD |3] and high-T c superconductivity We 
note that if we are able, here, to establish vorton stability in vacuum (i.e. 
resisting collapse due to the powerful tension of a relativistic string), then 
these results bode well for their as yet untested stability in less extreme phys- 
ical situations. We emphasise that here we are focussing on classical vorton 
stability, while for a discussion of quantum stability the reader is referred to 
[H] and references therein. 

The structure of this letter is as follows: After summarizing our previous 
results for vorton equilibrium states, we will examine additional relevant 
perturbations which we have observed as a result of our work. Finally, we 
will give an account of our numerical results on vortons and their long-term 
stability given elliptically perturbed initial conditions. 

2 The framework 

Our study is based on the model first proposed in ref. [2] : 

£ = {d^){d^) + + {d,a){d^Y - ^{\^ - V lf 

-^(w\ 2 -vi) 2 -m\ 2 w\\ a) 

where the two complex scalar fields 4> and a are minimally coupled, each 
being invariant under 2/(1) transformations. 

If the constants of the theory (the couplings and the vacuum expecta- 
tion values) are chosen carefully, one can break the 0-symmetry, leading to 
|0| = r]^, while keeping a = in the vacuum, because of the non-vanishing 
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interaction term. Along a cosmic string, the interaction vanishes, and a 
can form a condensate. (For a detailed acount on this model, see ref. [3]) 
We have been able to give necessary conditions on the parameters of the 
lagrangian for the string to be superconducting: 

A*<3§, vl< V j- ■ (2) 

This \a\ condensate can also carry charge and currents along the string 
(taken to lie on the 2-axis), as can be readily seen from the ansatz: 

a = \a\(r)e i{ujt+kz) , (3) 

which induces a (Noether) charge Q and a current J on the string worldsheet: 

Q = u [ dz [ dS H 2 , J = k [ dz [ dS \a\ 2 , (4) 



as well as a topologically conserved quantity, the winding number: 

X = Jdz±. (5) 

These features alter significantly the standard string cosmology scenario, 
with perhaps the most striking consequence of superconductivity being the 
classical stability of string loops. These loops, which were dubbed vortons 
in ref. [3 HI, cannot decay because of the angular momentum of their charge 
carriers. 

We present in figUJa typical, straight superconducting string profile. The 
parameters used here saturate the bounds given by eq.JSJ), in order to get a 
condensate well localised along the string, appropriate for numerical simu- 
lations with boundary restrictions. In general, one can expect a condensate 
somewhat broader than the underlying vortex: 5 a > 5^, or equivalently, 
Pvl ~ |A CT ?7 2 = m 2 < m 2 = ^X^r] 2 . It is probably worth emphasizing here 
the model dependence of vorton studies: the parameter space is very broad, 
and we need a specific particle physics model fixing the couplings to make 
more precise predictions about their properties. 

In our previous paper, we have been able to analyticaly characterize the 
equilibrium states proving, in principle, that vortons should occur for every 
initial non-zero value of Q and N. If, as usual, we call fi the string tension, 
and we define 

f dS\a\ 2 , S 4 = / dS\a\\ (6) 
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Figure 1: Profiles of and a, plotted in units of and r\ a \ the parameters 
here are = 1.5, rj^ = 1.0, \ a = 10.0, r] a = 0.5, (3 = 1.5 



we can recall the following equilibrium conditions. In the chiral case u 2 = k 2 , 
vortons will shrink or expand until: 



uj 2 = k 2 = fa - iA ff S 4o )/(2S ) 



(7) 



where the suffix o denotes the chiral value of a quantity (or equivalently, 
when uj — k — 0). In the electric {uj 2 > k 2 ) or magnetic (uj 2 < k 2 ) regimes, 
given that Eqjv = Q/N, we showed that the vorton state would minimize 
the following function of uj 2 — k 2 



u 2 : 



£ = N 



(A* 



QN 



E(u 2 ) 



Z 2 (u 2 ) 



(8) 



which clearly admits such a minimum. Given an initial state with small seed 
charges and currents, we demonstrated that the smaller final equilibrium 
state will generically be located away from the chiral state (which is not an 
attractor - see ref. |3] for details). 
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3 Stability analysis 



We employed numerical methods to investigate the behaviour of these equilib- 
rium vorton states, which enabled us to identify and analytically characterise 
various instabilities corresponding to different regimes of the loop. We have 
already proved in [3] that a straight superconducting string cannot evolve too 
deeply into the magnetic regime (u 2 < k 2 ), because of nonlinear effects that 
force the condensate to become pinched locally and to unravel, thus losing 
winding number. 

Our analytical analysis, and a suitable ansatz for £ and E 4 , gave a precise 
criterion for the threshold value ki nst of this instability. Let's define a = 
^S 4o /S , and k c , the maximum k value associated with a non-vanishing 
condensate. Then, 

111 
VT~ - ~ + VI ( 9 ) 

K inst a ft c 

We have also studied the chiral and electric cases, and we were unable 
to find an instability of this kind, suggesting that these regimes are stable 
against a 'pinching' perturbation. There is, however, another potential in- 
stability particularly relevant in the electric regime. 

In the model presented in eq. (JTJ, the radial loop equilibrium discussed 
above in [HI is only valid if one assumes that the charge and the current remain 
localised on the string's worldsheet. When we consider small loops, especially 
those tractable numerically, this assumption is not necessarily valid. The 
potential cost of the vortex/condensate separation can be evaluated using the 
usual profiles for the vortex, and we obtain the following stability criterion 
(see ref. 

where R is the radius of the vortex loop. Generically, the right hand side of 
(fTtij) is quite large, and so stable loops require R » 5 a . This is a problem 
for real vortons in the electric regime. The width of the condensate diverges 
as they go deeper into the electric regime, and so these loops can prove quite 
hard to stabilise. Numerically, memory limitations are such that we cannot 
achieve more than R ~ 205 a , and thus we are faced with this potential 
splitting instability even in the chiral regime, since this is just the order of 
magnitude given by (fT0|) . as we have demonstrated in our simulations. 

To overcome this difficulty, we have modified the interaction term in 
eq.JH), Vi n t = /3|0| 2 |ct| 2 , by considering a toy model with V int = /3'|0| 6 |cr| 2 . 



f*S§9> (10) 
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Figure 2: Profiles of and a in the modified Witten model, plotted in units of 
r]^ and r] a ; the parameters here are = 0.5, t]^ = 1.0, A CT = 18.0, t] a = 0.35, 
with the new interaction parameter (3' = 3.3 

Because of the higher power of involved, the effective potential seen by 
the condensate is much closer to a square well, which allows a to build up 
a higher and more robust condensate, as can be seen from figEl This repre- 
sents only a quantitative, rather than qualitative, modification of the vorton 
model with the new criterion for the splitting instability becoming less strin- 
gent. Then, the constraint given by (fT0|) is relaxed to R/5 a < 5, and we can 
create numerical configurations in which the condensate and the underlying 
vortex are held together tightly. 

4 Vorton simulations 

The theory defined in the previous section is ideal for numerical simulations. 
Our code evolves the fields according to the full 3D equations of motion aris- 
ing from eq.fJH), using a lattice-inspired hamiltonian formalism [H]] modified 
from |11| . energy and charge conservation in all simulations was maintained 
to below 1% accuracy. 

First, we study the perfectly circular chiral case. The initial configuration 
was obtained using an SOR relaxation routine to calculate the radial profiles. 
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t = o t = 1000 



Figure 3: Vorton simulation, showing \(f>\ (above) and the real part of cr, at 
t = and t = 1000 (note the uniform winding of cr, meaning a uniform 
current) . 

We then used the ansatz: 

cr(r,6A) = W\(r,9)e^ + ^, (11) 

which describes a homogeneous chiral vorton. Here, we have neglected small 
corrections due to the curvature. 

We then let the loop evolve with our code, using (Dirichlet) reflective 
boundary conditions which do not act to stabilise the configuration. We 
could observe the loop slowly oscillating around its equilibrium position, in 
agreement with the radial analysis we have given (see also [121 EH EH CL5J). 
As can be seen in figGl the whole structure appears to be remarkably stable. 

We now turn our attention to vorton stability with respect to perturba- 
tions in the eccentricity e, by considering a loop with e < 1. To ensure that 
the current is initially homogeneous, we have to consider the modified ansatz: 

a(r,9,s) = \a\(r,9)e ks+tJt . (12) 

where s is the arclength along the string. This arclength is given by an 
elliptic integral which we evaluate accurately using the Gauss- Tschebychev 
algorithm (as can be found in e.g. |TK]). 
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Our observations are summarized in fig. 0J where we can see clearly that 
the loop is actually oscillating between its initial configuration, and another, 
somewhat larger loop with higher eccentricity, tilted in the direction opposite 
to the current flow. We note also that this behaviour satisfies T ~ L, where 
T is the period of the movement, and L is the length of the loop. 

This can be understood if one considers the effect on the current of squash- 
ing a circular vorton, to make it look like an ellipse. To lowest order, the net 
effect on the angular momentum of the current is to induce a correction of 
the form: 

M = M (l+ e -). (13) 

Thus, the angular momentum of the current is incresaed with the eccentricity, 
and this has to be compensated by a rotation of the loop in the opposite 
direction (we discuss the transfer of momentum at length elsewhere |9j). 
Despite the oscillating eccentricity, ultimately these loops tend to evolve 
towards more circular configurations, as can be seen from the last two plots 
in figH 

Eccentric loop configurations retained their identity for more than 10,000 
time steps (many light-crossing times), and so these simulations appear to es- 
tablish that stable vortons should form during the evolution of the universe. 
Two small caveats to this conclusion remain. First, in the very longest simu- 
lations of over 30,000 timesteps, the build-up of background radiation (due to 
the reflective boundary conditions) causes some friction on the time-varying 
current which is eventually driven towards the less stable magnetic regime. 
We are developing absorbing boundary conditions for massive radiation to 
test the significance of this boundary artifact. Secondly, there is the sensitiv- 
ity of these objects to their initial conditions, since a slight mismatch in the 
phase of the condensate can have dramatic consequences on the subsequent 
evolution of the loop. This problem is still under active investigation and 
will be discussed at greater length elsewhere. 

5 Conclusions and discussion 

We have presented in this paper the first simulations of vortons, which appear 
to strongly indicate their stability, and hence their cosmological relevance (as 
well as in other physical contexts). Stable vortons may have profound con- 
sequences for cosmology. These objects can be very massive, and they could 
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Figure 4: Plots of the zeros of in the transverse plane of an elliptical chiral 
vorton: the first four plots show the f?rst period, and the last two are taken 
after two and three periods (initially, e = 0.4). 



account for dark matter, or even dominate the universe [13 HE). They are 
naturally associated with very high energies, and this makes them prime can- 
didates for high-energy astrophysical puzzles, like cosmic rays [TH] or gamma 
ray bursts [20]. Further study is required to test these hypotheses, but there 
can be little doubt that a deeper understanding of the microphysics of su- 
perconducting strings will further the confrontation between vortons and 
observations. 
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